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Abstract 

In the description of the production properties of gauge bosons (W ± , Z°, 7*) 
at colliders, the lowest-order graph normally is not sufficient. The contribu- 
tions of higher orders can be introduced either by an explicit order-by-order 
matrix-element calculation, by a resummation procedure or by a parton- 
shower algorithm. Each approach has its advantages and disadvantages. 
We here introduce a method that allows the parton-shower algorithm to 
be augmented by higher-order information, thereby offering an economical 
route to a description of all event properties. It is tested by comparing with 
the p± spectrum of W bosons at the Tevatron. 
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The W ± and Z° bosons have been extensively studied at colliders, in order to test the 
standard model [|l|. In recent years they have also made their debut as backgrounds to 
other processes of interest: top studies, Higgs searches, and so on. Here it is often the 
association of the W/Z with one or several jets that is the source of concern. Such higher- 
order corrections to the basic processes also serve as tests of QCD. It is therefore of some 
interest to improve the accuracy with which gauge boson production can be described. 

In this letter we will take the W ± production at hadron colliders as a test bed to develop 
some ideas in this direction. Specifically, we will discuss how to improve the lowest-order 
description qq' — > W ± by a merging of the first-order matrix elements qq' — > gW ± and 
qg — ► q'W 1 * 1 with a leading-log parton shower. However, the formalism is valid for all 
colourless massive vector gauge bosons within and beyond the standard model: 7*, Z°, 
Z'°, W /=t , and so on. It also applies e.g. in e + e~ — > 7Z . One could in addition imagine 
extensions to quite different processes, such as Higgs production by gg — > h°, but this 
would require further study. 

The outline of the letter is the following. First we discuss various approaches to W 
production, and their respective limitations. Then we zoom in on the shower method and 
introduce a matrix-element-motivated method to improve it. Finally we compare with 
data, specifically the W transverse-momentum spectrum at the Tevatron, and draw some 
conclusions. 

In essence, one may distinguish three alternative descriptions of W production: 

1. Order-by-order matrix elements. By a systematic expansion in powers of a s , a quite 
powerful machinery is obtained. For instance, there are calculations of the total 
Drell-Yan cross section to second order and for the associated production of a 
W and up to four partons at the Born level ||. A main problem is that there is no 
smooth transition between the different event classes, as one parton becomes soft or 
two partons collinear. The method is therefore better suited for exclusive questions 
than for an inclusive view of W event properties. 

2. Resummed matrix elements. Here the effects of multiple parton emission are re- 
summed, in impact-parameter or transverse-momentum space |4[]. Inclusive quan- 
tities such as the p±w spectrum can be well described in such an approach, but it 
should be noted that a nonperturbative input is required. The standard formal- 
ism does not give the exclusive set of partons accompanying the W, however, and 
internally does not respect correct kinematics. 

3. Parton showers. The parton-shower approach generates complete events, with cor- 
rect kinematics. An arbitrary number of partons is obtained, with a smooth and 
physical transition between event classes ensured by the use of Sudakov form factors. 
On the other hand, the shower approach is formally only to leading-log accuracy 
(although many detailed choices are made to maximize agreement with next-to- 
leading- log results), and the description of the rate of exclusive parton configurations 
may be poor. 

Finally, note that the perturbative partonic stage is not observable in experiment, but 
instead the hadronic jet one. Traditional hadronization descriptions, such as string frag- 
mentation 0, are intended to be universal if applied at some low cut-off scale Q ~ 1 GeV 
of the perturbative phase. This perfectly matches the shower approach, but causes prob- 
lems in the use of matrix elements. 

Given their complementary strengths, it is natural to attempt a marriage of the matrix- 
element and parton-shower methods, where the rate of well-separated jets is consistent 
with the former while the substructure of jets is described by the latter. The simpler 
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solution, matching, is to introduce a transition from one method to the other at some 
intermediate scale || [7], §]. Such an approach is convenient for descriptions of exclusive 
jet topologies, but tend to suffer from discontinuities between event classes and around 
the transition scale. More ambitious is the merging strategy, where matrix-element in- 
formation is integrated into the shower in such a way as to obtain a uniform and smooth 
description. This approach so far has only been implemented for the merging to 0(a s ) 
of e + e~ — ► qq with e + e~ — > qqg || 0. We will here introduce a corresponding 0(a s ) 
merging in hadronic W production. Further details may be found in JIO . 



Since we neglect the decay of the W, alternatively imagine it decaying leptonically, all 
QCD radiation occurs in the initial state. We will base our approach on the initial-state 



shower algorithm of [ITJ, as implemented in Pythia [fL2| . The principle of backwards 
evolution implies that a shower may be reconstructed by starting at the large Q 2 scale of 
the hard process and then gradually considering emissions at lower and lower virtualities, 
i.e. earlier and earlier in the cascade chain (and in time). 



The starting point is the standard DGLAP evolution equation [13 



-W 1 !* ^ /„(*',*) p^m, (i) 



dt a Jx x' 2n 

with fi the distribution function of parton species i, x the momentum fraction carried by 
the parton, t = ln(Q 2 /AQ CD ) the resolution scale, and P a ^i, c {z) the AP splitting kernels 
for parton b obtaining a fraction z = x/x' of the a momentum. Normally the evolution 
is in terms of increasing t, but in the backwards evolution t is instead decreasing. Then 
the DGLAP equation expresses the rate at which partons b of momentum fraction x are 
'unresolved' into partons a of fraction x', in a step dt backwards. The corresponding 
relative probability is dPb/dt = (1//&) (dfo/dt). The probability that b remains resolved 
from some initial scale t max down tot < t max is thereby obtained by a Sudakov form factor 

S b {x,t;t max ) = exp - / — — — — — dt 



t f b (x,t>) dt' 

I /" w a* f 1 dx ' a ^ jgfcVO p t \\ 

= exp [- I dt £ ^ - — P a ^ c{ z) J 

= exp [- 1 dt ? L dz -^umT p - Uz) ) ■ (2) 

From this expression it is a matter of standard Monte Carlo techniques to gener- 



ate the complete branching a — > be e.g., the t distribution of the branching is 
— dSb(x, t;tmax)/dt Given parton a, one may in turn reconstruct which parton branched 
into it, and so on, down to the starting scale Qq. In each branching, the t scale gives the 
tmax value of the branching to be considered next, i.e. the Q 2 values are assumed strictly 
ordered. 

The definition of the Q 2 and z variables is not unambiguous. Referring to the notation 



of Fig. [IJ and to the branching 3^1+4, the Q 2 scale in our algorithm [IT] is associated 
with the spacelike virtuality of the produced parton 1, Q 2 = —p 2 , while z is given by the 
reduction of squared invariant mass of the contained subsystem, z = (pi + P2) 2 / (j>3 + P2) 2 ■ 
In the limit of collinear kinematics, Q 2 = 0, one recovers the momentum fraction z = 
Pi/p3- The z definition couples the two sides of the events, so that the order in which the 
branchings 3 — > 1 + 4 and 5^2 + 6 are considered makes some difference for the final 
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Figure 1: Schematic picture of an initial-state parton shower, extending from both sides 
of the event in to the W. 



configuration. The rule adopted is therefore to reconstruct branching kinematics strictly 
in order of decreasing Q 2 , i.e. interleaving emissions on the two sides of the event. 

Now let us compare the step from qq' — > W to qq' — > gW between the matrix- 
element and parton-shower languages. Since only one branching is to be considered, the 
comparison has to be with a truncated shower, e.g. where only the branching 3^1+4 
occurs in Fig. [1]. The 2^2 process thus is q(3) + q'(2) — > g(4) + W(0), for which 
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The matrix element for qq' — > gW can be written as [14 



da 
~dl 



ME 



§ 2tt 3 



u 2 + 2m^s 



tii 



(4) 



Here ao is the cross section for qq' — ► W, a = (ii 2 a era /3 sm 2 9w^ 2 w)\V q q , \ 2 S(l—rn 2 w /xiX2s) 
in the narrow-width limit, with 5(1 — m 2 N /xiX 2 s) t— > J dz5(l — m 2 ^/zx 3 X2s) in the 2^2 
process kinematics. (The details of the <7o factor are not relevant for the point we want 
to make, so the presentation is intentionally sketchy.) Now rewrite eq. (^) in terms of z 



and Q , using eq. 
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(5) 



PSl 



We here easily recognize the splitting kernel for q — > qg, i.e. the matrix element reduces 
to the the normal shower expression in the collinear limit, as it should be. Some extra 
but trivial work is necessary to include the convolution with parton distributions, which 
involves fi(xi,Q 2 ) in lowest order and f-z(x 3 ,Q 2 ) for the 0(a s ) processes. 

In order to study how the shower populates the phase space, it is straightforward to 
translate back the above expression, 



da 
~dt 



a a s 4 s 2 



m 



w 



PSl 



s 2n 3 t(t + u) 



(6) 



To this we should add the other possible shower history, where the gluon is emitted 
by a branching 5—^2 + 6 instead; after all, the matrix-element expression contains both 
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amplitudes. The collinear singularity Q 2 — > here corresponds to emission along direction 
2 rather than direction 1. In that case the roles of t and u are interchanged, and the cross 
section da/dt|ps2 is easily obtained. The total shower rate is given by the sum, 



da 
dt 



PS 



da 
dt 
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PSl 
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s 2n 3 tw 



(7) 



Thus the singularity structure of the parton-shower and matrix-element rates agree, giving 
a ratio 

(da/dt) M E P + M 2 + 2m^ 2tu 



-Rqq'->gw(s, t) 



constrained to the range 
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The same exercise may be carried out for qg — > q'W: 
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Note that, unlike the qq' — > gW process, there is no addition of two shower histories when 
comparing with matrix elements, since here also the latter contains two separate terms 
corresponding to qg and gq initial states, respectively. 

The qq' — > gW process receives contributions from two Feynman graphs, t-channel 
and li-channel, and the shower thus exactly matches this set, although obviously it does 
not include interference between the two. The qg — > q'W process is different, since only 
its u-channel graph is covered by the parton-shower formalism, while the s-channel one 
has no correspondence. Since this latter graph is free from collinear singularities, the 
shower is not misbehaving in any regions of phase space because of this omission, but it 
is interesting to speculate that the larger value for i? qg ^ q / W (s, t) than for -R q q'_ tg w('5, t) 
partly may have its origin here (remember that a larger R(s, t) means a smaller shower 
emission rate). 

Based on the above exercise, the standard parton-shower approach may be improved 
in two steps. The first is to note that, since the shower so closely agrees with the correct 
matrix-element expression — much better than one might have had reason to expect - 
it is safe to apply the shower to all of phase space, i.e. to have QiL v ~ s rather than 
the more traditional shower-generator limit Q 2 ^ 
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The older choice was 

2 ^ ™2 



inspired in part by the fear of a completely erroneous behaviour for Q 3> ?r%, in part by 
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the typical factorization scale used for parton distributions in W cross-section formulae. 
Such a scale choice can be motivated by doublecounting arguments. Most easily this is 
seen in pure QCD processes, where a 2 — > 3 process such as gg — ► qqg could be obtained 
starting from several different 2 — > 2 processes, and classification by the hardest (most 
virtual) subgraph is necessary to resolve ambiguities. Correspondingly, a W production 
graph could be reclassified once some parton has Q 2 > m 2 ^. That is, in general, one 
would have to consider QCD processes where the emission of a W is allowed as a 'parton 
shower' correction. (The s-channel graph in qg — > q'W is an example of this kind.) 
Doublecounting is not an issue, however, once we decide to represent the full W cross 
section by qq — > W. (Remember that the shower does not change total cross sections.) 

The second step is to use standard Monte Carlo techniques to correct branchings 
in the shower by the relevant ratio R(s,t), to bring the shower parton-emission rate in 
better agreement with the matrix-element one. This correction is applied to the branching 
closest to the hard scattering, i.e. with largest virtuality, on both sides of the event, i.e. 
for 3^1 + 4 and 5 — > 2 + 6 in Fig. |l| By analogy with results for time-like showers 
0, one could attempt to formulate more precise rules for when to apply corrections, 
but this one should come close enough and is technically the simplest solution. (For 
instance, while our cascade is ordered in Q 2 rather than in p\, the emission with largest 
p 2 ± normally coincides with the largest Q 2 one, so either criterion for when to apply a 
correction would give very similar results.) For a q — > qg shower branching, where the 
correction factor R(s,t) = R q q'^gw(S,t) < 1, a candidate branching selected according 
to the Sudakov factor in eq. (fj) is then accepted with a probability R(s,t). In case of 
failure, the evolution downwards in Q 2 is continued from the scale that failed (the 'veto 
algorithm', ensuring the correct form of the Sudakov). For a g — > qq branching, the fact 
that R{s,t) = -R q g_ +q 'w(S, t) > 1 means that the procedure above cannot be used directly. 
Instead the normal g — > qq branching rate is enhanced by an ad hoc factor of 3, and the 
acceptance rate instead given by R(s,t)/3 < 1. 

Even with this injection of matrix-element information into the parton shower, it is 
important to recognize that the shower still is different. The hardest emission is given 
by the matrix-element expression times the related Sudakov form factor, thus ensuring a 
smooth p± spectrum that vanishes in the limit p± — ► 0. By the continued shower history 
(without any matrix-element corrections), further emissions pick up where the Sudakov 
factor suppresses the hardest one, giving a total p± spectrum of emitted partons that is 
peaked at the lower p± cut. This total spectrum is similar to the matrix-element one, 
but deviates from it in that the shower includes kinematical and dynamical effects of 
gradually having partons at larger and larger x values and possibly of different species 
at each softer emission. In some respects, it thus provides a more sophisticated approach 
to resummation for the properties of the recoiling W. It also gives exclusive final states, 
including the possibility for the emitted partons (such as 4 and 6 in Fig. |TJ) to branch in 
their turn. 

The parton shower redistributes the W's in phase space but does not change the total 
W cross section. It is thus feasible to use a higher-order calculation of this cross section 
as starting point, although we did not do it here. If higher orders enhance the total cross 
section by a factor K (with K a function e.g. of rapidity) relative to the lowest-order 
qq' — ► W one, the implication of eqs. (H) and fllCf ) is that the one-jet rate is enhanced 
by the same factor. If instead another factor K' is wanted here, the respective R(s, t) 
weight could then be modified by a factor K' / K. Note, however, that it is more difficult 
to introduce such a K'/K factor consistently, since it presupposes a common definition 
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Figure 2: The p^w distribution in pp collisions at 1.8 TeV. Parton distributions and 
a s are frozen and only one emission at a time is allowed in the shower. Events are 
classified either as (a) qq' — > gW or (b) qg — > q'W. 



between showers and matrix elements of what constitutes a jet. 

As a first check, we want to confirm that the shower algorithm works as intended, 
reproducing the matrix-element expressions where it should. Thus the shower is artificially 
modified so as only to generate one branching at a time. In order to eliminate the influence 
of the Sudakov and the change of kinematics by previously considered emissions, the 
shower is restarted from each Q 2 actually selected above the cut, but returning to the 
original kinematics for qq' — > W. Furthermore parton distributions and a s are frozen, so as 
to avoid any scale-choice mismatches. The resulting W transverse momentum spectrum 
is shown in Fig. ^, classified by the two possible branchings. In the old scheme, with 
Qmax = m W' the drop of the p_i_w spectrum at p_i_w ~ m w is easily visible. Already 
the modification to Q^ax = s 'intermediate' curves) brings a marked improvement, 
and the further introduction of the R(s, t) weighting ('new') results in good agreement 
between the shower and the matrix elements. 

The R(s, t) factors are further studied in Fig. |3|. It is seen that R(s, t) is close to unity 
for most of the branchings (note the logarithmic scale). Also that R(s, t) — > 1 as p± — > 0, 
in accordance with the demonstrated agreement of the parton shower and matrix elements 
in the collinear limit. At large p± values, the R(s, t) factors enhance the importance of 
the qg — ► q'W process relative to the qq' — ► gW one by about a factor of 2. When the 
two processes are not separated, the partial cancellation of having one R(s, t) a bit above 
unity and the other a bit below leads to a rather modest net correction to p± spectra. 

We now present results for the new parton shower, in all its complexity, i.e. with 
normal complete showers augmented by the two-step correction process described above. 
In Fig. |] the shower p_i_w distribution is compared with experimental data from the 
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Figure 3: R(s,t) distributions in pp collisions at 1.8 TeV. (a) The inclusive distri- 
bution, (b) The average value as a function of pj_w- Events are classified either as 
qq' — > gW or qg — > q'W. 



DO collaboration JyJ. The agreement is good for large p±w, but the shape at small 
p_i_w is rather different, with less activity in the shower than in data. These results are 
for the default Gaussian primordial k±_ spectrum of width 0.44 GeV, as would be the 
order expected from a purely nonperturbative source related to confinement inside the 
incoming hadrons. By now several indications have accumulated that a larger width is 
needed, however |]16f , although the origin of such an excess is not at all understood. 
One hypothesis is that some radiation is overlooked by an imperfect modelling of the 
perturbative QCD radiation around or below the Qo cut-off scale. Whatever the reason, 
we may quantify the disagreement by artificially increasing the primordial k±_ width. 
Fig. [| shows that an excellent agreement can be obtained, at all p_i_\y values, with a 
4 GeV width. The p±~w distribution is essentially unchanged at large values, i.e. only 
the region p^w ~ 20 GeV is affected. In order to put the 4 GeV number in perspective, 
it should be noted that this is introduced as a 'true' primordial k±, i.e. carried by the 
parton on each incoming hadron side that initiates the initial-state shower at the Qo scale. 
If such a parton has an original momentum fraction xq and the parton at the end of the 
cascade (that actually produces the W) has fraction x, the W only receives a primordial 
k± kick scaled down by a factor x/xq- In the current case, and also including the fact 
that two sides contribute, this translates into a rms width of 2.1 GeV for the primordial 
k± kick given to the W. This number actually is not so dissimilar from values typically 
used in resummation descriptions 

In summary, we see that it is possible to obtain a good description of the complete p±w 
spectrum by fairly straightforward improvements of a normal parton-shower approach. 
Corresponding improvements can also be expected for the production of jets in association 
with the W. Especially, the good matching offered to hadronization descriptions in this 
approach allows a complete simulation of the final state, including the addition of a 
(possibly process-dependent) underlying event. This way it should be possible to address 
e.g. the ratio of events with/without a jet accompanying the W, where CDF and DO have 
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Figure 4: Transverse momentum spectrum of the W; full parton-shower results 
compared with data from the DO Collaboration ]T5| . (The error bars include both 
statistical and systematic uncertainties.) 



obtained partly conflicting results, however using different jet definitions [[□]]. 

We end by reiterating that the formalism presented here is universal, in the sense 
that the formulae in this paper are not unique for the W, but shared by all vector gauge 
bosons, after an appropriate replacement of m-yv and the constants in the gq prefactor. 
Specifically, the reweighting factors R(s, t) need only be modified to reflect the mass of the 
current resonance. The method therefore should offer an accurate and economical route 
to the prediction of kinematical distributions for a host of new particles, to be searched 
for at the Tevatron and the LHC It is also likely that similar approaches can be developed 
for other classes of processes. 
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